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Abstract 

In this paper, the equations governing the unsteady flow of a perfect 
polytropic gas in three space dimensions are considered. The basic sim- 
ilarity reductions for this system are performed. Reduced equations and 
exact solutions associated with the symmetries are obtained. This results 
is used to give the invariance of system up to Galilean motions of space- 
times K 4 . Then, an optimal system of one-dimensional sub-algebras for 
symmetry algebra of this system is given. 
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1 Introduction 

The equations governing the unsteady flow of a perfect polytropic gas in 
three space dimensions are 

du du du du 1 dp 

dt dx dy dz q dx 

dv dv ^dv ^dv 1 dp _ ^ 

dt dx dy dz q dy 

dw dw dw dw 1 dp , . 

dt dx dy dz q dz ' 

dq ( du dv dw \ dq dq ^ dq _ ^ 

dt ^\dx dy dz ) dx dy dz 

dp / du dv dw \ dp ^ ^ dp ^ dp _ ^ 

dt ^ v dx dy dz ) dx dy dz 

where t is the time and x, y, z are the space coordinates; q(t, x, y, z) is the 
density, p(t,x,y,z) the pressure, u(t,x,y,z), v(t,x,y,z) and w(t,x,y, z) 
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the velocity components in the x, y and z directions, respectively, and 7 £ 
M is the adiabatic index. These equations are called "three dimensional 
Euler equations of gas dynamics" [8] and section 6.3 of [9] . 

In this paper, we consider the equations governing the unsteady flow 
of a perfect polytropic gas in three space dimensions. The basic similarity 
reductions for this system are performed. Reduced equations and exact 
solutions associated with the symmetries are obtained. We find an opti- 
mal system of one-dimensional sub-algebras for symmetry algebra of this 
system. 

In section 2 we find the full symmetry algebra of system (1). Chapter 
3 is devoted to the group-invariant solutions to the system. The structure 
of full symmetry algebra illustrated in section. The last section deals with 
a optimal system of sub-algebras. 

2 Lie symmetry of the system 

The classical method for finding symmetry reductions of PDE is the Lie 
group method of infinitesimal transformations. To apply the classical 
method to (I), we consider the one-parameter Lie group of infinitesimal 
transformations in (t,x,y,z,u,v,p,q) given by 



t 


= t + s.&(t,x,y,z, 


u, 


v,p,q) + 0(s 2 ), 


X 


= x + s.&(t,x,y,z, 


u, 


v,p,q) + 0(s 2 ), 


V 


= y + s.£ 3 (t,x,y,z, 


u, 


v,p,q) + 0(s 2 ), 


z 


= z + s.£ 4 (t,x,y,z, 


u, 


v,p,q) +0(s 2 ), 


u 


— u + s.r/i (t, x, y, z. 


, u. 


.v,p,q) + 0(s 2 ), 


V 


= v + s.rj 2 (t,x,y,z. 


u, 


v,p,q) + 0(s 2 ), 


w 


= w + s.r/ s (t,x,y,z 


,u 


,v,p,q) + 0(s 2 ), 


p 


= p + s.f] 4 (t,x,y,z, 


u, 


v,p,q) + 0(s 2 ), 


q 


= q + s.r) 5 (t,x,y,z, 


u, 


v,p,q) + 0(s 2 ). 



where s is the group parameter. One requires that this transformation 
leaves the sub-manifold S invariant. Let 

X = £1 d t + £2 d x + £3 d y + £4 d z + ??i d u + ??2 d v + m 9 W +V4.d p + r) 5 d q 

be the corresponding infinitesimal transformation; where coefficients are 
real-valued C°° —functions of t, x, y, z,u,v, w, p, q. The infinitesimal crite- 
rion of invariance (1) will thus involve t,x,y, z, u, v, w, p, q and the deriva- 
tives of u,v,w,p,q with respect to t,x,y,z, as well as i = 1,2,3,4, and 
Tjj, j = 1, 2, 3, 4, 5, and their partial derivatives. After eliminating any de- 
pendencies among the derivatives of the m's caused by the system itself, 
we can then equate the coefficients of the remaining unconstrained partial 
derivatives of u to zero. This will result in a large number of elementary 
partial differential equations for the coefficient functions £», r/j of the in- 
finitesimal generator, called the determining equations for the symmetry 
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group of the given system. In this case, we find the determining equations 
for the symmetry group of the system (1) to be the following: 

£i,t 2 = £i,x = £1,1/ = £i,z = = £i,« = £1,10 = £1,9 = £i, P = 0, 

£2,2; = £,2,u = &,v = £2,™ = £2,9 = £,2,p = £ 2 ,t 2 = £2,ty = £2, y 2 = 0, 

^3, y = £3,11 = £3,1) = £3,™ = £3,9 = £3,p = C3,t 2 = 0, (3) 

£4,2 = £,A,w = £,A,q = £4,p = 0, 

?4,t 2 = £,4,,tx = £4,43/ = £4, x 2 = £,4,xy = t,4,y 2 = 0, 

£,3,x = £2, y, £3,z + £4, a = 0, £4,1, = £4,11, £2, z + £4,1 = 0, 

and 

Til = -«£l,t + £2,t + <2,y - w£4,^, 
»J2 = ~W^4,y + £3,t - «£2,s, - «£l,i, 

??3 = «£4,x +£4,t + v(,A,y ~ «>£l,t, (4) 

»74 = qrfc/p + 2q^ t , 

V5,t = V5,x = r?5, H = V5,z = r?5,u = V5,v = »75,«i = »75,g = 0, 
T)5 = pri5, P - 

First, equations (3) require that £iS are just 
£1 = ait + a,2, 

£2 = a 3 t + a t x + a 5 y + a 6 z + a 7 , 

£3 = a$t — a^x + ci4j/ + a 9Z + aio, (5) 
£4 = ant — ciax — agy + cuz + ai2, 

where ai,---,oi 2 are arbitrary constants. Then, equations (4) and (5) 
requires that 

771 = (04 - ai)u + a 5 t> — a 6 w + a 3 , 

r/2 = — a&u + («4 — ai)v + (14W + as, 

r]3 = —a^u — agv + (04 — ai)w + an, (6) 

f]4 = (2oi - 2a 4 + 013)0, 

rfs = ai3P, 

where ai3 is an arbitrary constant. Therefore 
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Theorem 1. The Lie algebra of infinitesimal symmetries of the system 
(1) is spanned by the 13 vector fields 



(shifts of origin) 



x l 


= d x , 


x 2 


= 9y, 


x 3 


= d z , 


Xi 


= d t , 



Xs = t dx + d u , 

X§ — td y + d v , (uniformly frame motions) 

X 7 = td z + d w , 

(7) 

X s = t d t - u d u - v d v - w d w + 2q d q , 

X 9 = xd x + yd y + zd z + ud u + vd v + w d w — 2r d r , 



Xio = y d x - x d y + v d u - u d v , 

Xn = — zdy +yd z — w d v +vd w , (rotations of reference frame) 
Xi2 = — zd x + x d z — w d u + ud w , 



Xi 3 = qd q + pd p , 

These infinitesimal symmetries will generate a Lie algebra Q over the field 
of real or complex numbers. The commutator table of Lie algebra g for (1) 
is given below, where the entry in the ith row and jth column is defined 
as [Xi,Xj] = Xi.Xj - Xj.Xi, i,j = 1, ■■■,13. 





Xi 


x 2 


x 3 


Xi 


x 5 


x 6 


x 7 


Xs 


x 9 


Xio 


Xn 


Xi 2 


Xl 3 


x 1 




























-x 2 





x 3 





x 2 


























x 2 


Xi 


X 3 








x 3 


























x 3 





-x 2 


-Xi 





x 4 














X! 


x 2 


X 3 


Xi 

















x 5 






















-X 5 X 5 


-x 6 





x 7 





x 6 











-x 2 











-x 6 


X G 


X 5 


x 7 








X 7 











-X 3 
















-X 6 


-X 5 





X s 











-Xi 


x 5 


X 6 


x 7 




















x 9 


-x 1 


-x 2 


-x 3 





-x 5 


-X 6 


-x 7 




















Xio 


x 2 










X 6 


-X 5 














— Xl2 


Xn 





Xu 





-X 3 


x 2 








-x 7 


x 6 








Xi 2 





—Xio 





X\2 


-x 3 





Xi 





-x 7 





x 5 








-Xn 


Xio 








X13 










































3 Reduction of the system 

Now, we look at the group- invariant solutions to the system (1). 

The system (1) is a sub-manifold M of jet space J 2 (R 4 ,R 5 ) defined 
by the following equations 

u t + uu x + vuy + wu z + p x /q = 0, 
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vt + uv x + vv y + wv z +p y /q = 0, 

w t + uw x + vw y + ww z + p z /q = 0, (8) 
qt+q (u x + v y + w z ) + uq x + vq y + wq z = 0, 
p t + jp (u x + v y + w z ) + up x + vp y + wp z = 0. 

Doing as section 3.1 of [6] and find (in a sense) the most general group- 
invariant solutions to the system (1). 

Theorem 2. The one-parameter groups gi(s) : M — > M generated by 
the Xi, i = 1, • • • , 13 are given in the following table: 

91 ■ (t,x + s,y,z,u,v,w,p,q) 

92 ■ (t,x,y + s,z,u,v,w,p,q) 
9s ■ (t,x,y,z + s,u,v,w,p,q) 
94 ■ (t + s,x,y,z,u,v,w,p,q) 

g 5 : (t,x + st,y,z,u + s,v,w,p,q) 
g 6 : (t,x,y + st,z,u,v + s,w,p,q) 

97 ■ (t,x,y,z + st,u,v,w + s,p,q) (9) 
gs : (e s t, x, y, z, e~ s u, e~ s v, e~ s w,p, e 2s q) 
g$ : (t,e s x, e s y,e s z, e s u,e s v, e s w,p,e 2s q) 
gw : (t, x coss + ysins, y cos s — a; sin s, z, 

«sins + u cos s,v coss — u sin s, w,p, q) 
gn : (t, x, y cos s — z sin s, z cos s + y sin s, 

u, v cos s — w sin s, w cos s + v sin s, p, q) 
gi2 : (t, x coss — z sin s, y, x sin s + z cos s, 

wcos s — w sins, w, w cos s + usin s,p, g) 
3i3 : (t,x,y,z,u,v,w,e s p,e s q) 

Where, entries give the transformed point 

exp(sXi)(t,x,y,z,u,v,w,p,q) = (t,x,y,z,u,v, w,p, q). 

Since each group gi, i = 1, • • • , 13, is a symmetry group of (1), then 

Theorem 3. Let i = 1, • • • , 13 ; s € K, and Ael - {0}. // 

u = U{t,x,y,z), v = V(t,x,y,z), w = W(t,x,y,z), 
p = P(t, x,y,z), q — Q(t, x, y, z), 
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is a solution of the system (1), so are the functions 

u (i) = u (i) {t,x, yi z), v (i) = V (l) (t,x,y,z), w w = W (l) (t,x,y,z), 
pW=pW{t,x,y,z), q (i) =QV(t,x,y,z), 



LUflCl fc. 






u« = U(t,x + s,y,z), 


= 


U(t,x,y + s,z), 


v (1) = V(t,x + s,y,z), 




V(t,x,y + s,z), 


W W = W(t,x + s,y,z), 


w m = 


W(t,x,y + s,z), 


= P(t,x + s,y,z), 


p (2) = 


P(t,x,y + s,z), 


— Oii t -1- q ii y\ 

q — d, -f b, y, zj . 


„(2) _ 
1 — 




u (3) = U(t,x,y,z + s), 


ti< 4 > = 


U(* + s,a;,y,z), 


= v(t,x,y,z + s), 


V W = 




w {3 > = W(t,x,y,z + s), 


(A \ 




p (3) = P(t,x,y,z + s), 


p (4) = 




g (3) = Q(t,x,y,z + s). 


,(4) = 


Q(t + s,x,y,z). 


u (5) = U(t,x + st,y,z) - s, 


u< 6 > = 


U(t,x,y + st,z), 


« (5) = V(t,x + st,y,z), 


„(6) = 


V(t, x, y + st, 2) — 


w (5) = W(t,x + st,y,z), 


w (6) = 


: W(t,x,y + st,z), 


p< 5 > = P(t,x + st,y,z), 


p (6) = 


P(t,x,y + st,z), 


q (5) = Q(t,x + st,y,z). 


g (6) = 


Q(t,x,y + st,z). 


u m = U(t,x,y,z + st), 




= X.U(X.t,x,y,z), 


= V(t,x,y,z + st), 


„(8) 


= X.V(X.t,x,y,z), 


w m = W(t,x,y,z + st) - s, 


w w 


= X.W(X.t,x,y,z) 


p {7) = P(t,x,y,z + st), 


p (8) 


= P^X.t, x, y, 2) , 


q (7) = Q(t,x,y,z + st). 


? (8) 


= jz-Q(X-t,x,y,z) 



tt (9) = ±.U(t,X.x,X.y,X.z), u (13) = U(t,x,y,z), 

v (9) = ±.v(t,A.a;,A.j/,A.z), « (13) = V (t, x, y, z) , 

w 0) = j-.w(t,X.x,\.y,\.z), w (13) = W(i,a;,y,«), 

p (9) = P(t,X.x,X.y,X.z), ' p (13) = X.P(t,x,y,z), 

= X 2 .Q(t,X.x,X.y,X.z). g< 13 > = X.Q(t, x,y, 2) . 



u ' — cos s.U[t, x cos s + y sins, y cos s — resins, zj 
— sin s.y (t, x cos s + y sin s, j/ cos s — a: sin s, 2^ 

i/ 10 ' = sin s.ll(t, x cos s + j/sin s, j/coss — :rsin s, 2) 
+ cos s.V(t, a; cos s + 2ysins,ycoss — a; sins, 2 



G 



,(10) 



,(10) 



.(10) 



W(t, x cos s + y sin s, y cos s — x sin s, 2) , 
P(t, x cos s + y sin s, j/ cos s — x sin s, 2) , 
Q(t, 2: cos s + j/ sin s, j/coss — £ sin s, 2) . 



(ii) 
(ii) 



„(") 



(11) 
.(11) 



U(t, x, j/cos s — z sins, 2 cos s + j/sin s), 
cos s.V(t, x, j/cos s — 2 sins, 2 cos s + j/sins) 
+ sin s.lWt, x, j/cos s — 2 sins, z cos s + j/sin s) , 
cos s.W(i, x, j/cos s — 2 sins, z cos s + j/sin s) 
— sin s.V (t, x, y cos s — 2 sin s, 2 cos s + j/ sin s) , 
P(t, x, y cos s — 2 sin s, 2 cos s + j/sin s) , 
Q(t,x,ycoss — 2 sins, 2 cos s + j/sins). 



,(12) 



,(12) 
,,(12) 



,(12) 
.(12) 



cos s.U (t, x cos s — 2 sin s, y, x sin s + 2 cos s) 
+ sin s.lWt, x cos s — 2 sin s, j/, x sin s + 2 cos s), 
V(i,xcoss — 2 sin s, j/, x sin s + 2 cos s), 
coss.lWt, xcoss — 2 sins, j/,xsins + 2 cos s) 
— sin s.U (t, x cos s — 2 sin s, j/, x sin s + 2 cos s) , 
P(i,xcoss — 2 sins, j/,xsins + 2COSs), 
Q(i, xcoss — 2 sins, j/,xsins + 2coss). 



4 Structure of Lie algebra g 

In this section, we determine the structure of full symmetry algebra Q. of 
system (1). 

The center 3 of g is Span R {Xi3}. Therefore, the quotient algebra 
fli ~ 2/3 is Span R {Yi, ■ ■ ■ , Y12}; where Y t := Xi + 3, i = 1, ■ ■ ■ , 12. 
The commutator table of Lie algebra Qi is given below, where the entry 
in the ith row and jth column is denned as = Yi.Yj — Yj.Yi, 
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i,j = !,-■■, 12. 





Fi 


Y 2 


Y 3 


Y 4 


Y 5 


Y 6 


Y 7 


Y s Y 9 


Y10 


Yn 


Yl 2 

























Yi 


-y 2 





Y 3 


y 2 























Y 2 


Yi 


Ys 





y 3 























y 3 





-Y 2 


-Yi 


y 4 














Yi 


Y 2 


Y 3 


y 4 











y 5 











-Yi 











-Y B Ys 


-Ye 





Y 7 


y 6 











-Y 2 











-y 6 y 6 


Ys 


Y 7 





Y 7 











-Ys 











-Yr Y 7 





-Y 6 


-Ys 


Y s 











-Yi 


Y 5 


Y 6 


Y 7 














Y 9 


-Yi 


-Y 2 


-Y 3 





-Y, 


-Ye 


-Y 7 














Y 10 


Y 2 


-Yl 








Y 6 


-Ys 











-Yi 2 


Y11 


Fn 





-Y 3 


Y 2 








-Y 7 


Y e 





Y12 





-Y10 


Y12 


-Y 3 





Yi 





-Y 7 





Y 5 





-Y11 


Y10 






The algebra gi is not solvable, because 
fli 1 ' :=[fli,fli] = Span R {Yi ) ...,Yio}, and < 2) := «] = fl« . 



Lie algebra Qi admits a Levi-decomposition into the semi-direct sum 

where r = Span R {Yi, • ■ ■ , Y9} is the radical of Qi (the largest solvable 
ideal contained in Qi) and S = Span R {Yio, Yn , Y12} is a semi-simple 
subalgebra of Qi. The radical r is uniquely defined but the semi-simple 
subalgebra S is not. S as an algebra is isomorphic to 3— dimensional special 
orthogonal algebra so(3) := {A € Mat (3 x 3;R) | A* = -A}, which is a 
simple algebra. 

The radical t is solvable, with and 

r (1) =Span R {Yi,...,Y 7 }, t (2) = Span R {Yi, Y 2 , Y 3 }, r (3) = {0}. 

It is semi-direct sum c = K Q 2 , where Q 2 = Span R {Yg, Yg} is isomor- 
phic to the two dimensional Abelian Lie algebra R 2 ; i.e. 

r = r (1) k R 2 . 

r' 1 ' is semi-direct sum = r (2) IX 03, where 03 = Span R {Y4, • • ■ , Y7} 
is isomorphic to the four dimensional Abelian Lie algebra R 4 , and r' 2 ^ is 
isomorphic to the three dimensional Abelian Lie algebra R 3 ; i.e. 

r (1) =R 3 k M 4 . 

Theorem 4. The symmetry algebra q of system (1) have the following 
structure 

g = R x (((M 3 k R 4 ) K R 2 ) k so(3)). 
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Theorem 5. There is a normal Lie-subalgebra of symmetry Lie-group 
of system (1) which is isomorphic to the Galilean group o/M 4 : 



Gal(4,l 




R G 0(3, M), s G R, and y, v G R 3 



Therefore, the system (1) is invariant up to Galilean motions of space- 
times R 4 (see [3J). 

Proof. Let b = Span R {Xi, ■ • • , X7, X10, Xn, X12}. Lie algebra structure 
of b is as ga[(3), the lie algebra of Galilean group Gal(4, R), and b an ideal 
of g. Therefore, there is a Lie-subgroup of G such that its Lie-algebra is 
b, by the Theorem 2.53 of Olver [6]. □ 

Conclusion 1. q is semi-direct sum of 

ga[(4) ^ Span R {Xi,---,X 7 ,Xi ,Xii,Xi 2 } 
and 3— dimensional Abelian algebra R 3 = SpanjXs, X9, X13}. 



5 Optimal system of sub-algebras 

As is well known, the Lie group theoretic method plays an important role 
in finding exact solutions and performing symmetry reductions of differen- 
tial equations. Since any linear combination of infinitesimal generators is 
also an infinitesimal generator, there are always infinitely many different 
symmetry subgroups for the differential equation. So, a mean of deter- 
mining which subgroups would give essentially different types of solutions 
is necessary and significant for a complete understanding of the invariant 
solutions. As any transformation in the full symmetry group maps a solu- 
tion to another solution, it is sufficient to find invariant solutions which are 
not related by transformations in the full symmetry group, this has led to 
the concept of an optimal system [7]. The problem of finding an optimal 
system of subgroups is equivalent to that of finding an optimal system of 
subalgebras. For one-dimensional subalgebras, this classification problem 
is essentially the same as the problem of classifying the orbits of the ad- 
joint representation. This problem is attacked by the naive approach of 
taking a general element in the Lie algebra and subjecting it to various 
adjoint transformations so as to simplify it as much as possible. The idea 
of using the adjoint representation to classify group-invariant solutions 
was due to [7] and [5] . 

The adjoint action is given by the Lie series 

Ad(exp(sX i )X j ) = Xj - spCi, Xj] + ^[X ; , [X;, Xj]] - ■ • • , (10) 

where [Xi,Xj] is the commutator for the Lie algebra, s is a parameter, 
and i,j = 1, • • • , 13. We can write the adjoint action for the Lie algebra 
Q, and show that 
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Theorem 6. A one-dimensional optimal system of (1) is given by 

1) X = X 8 + ^]i g a i X i , 

2) X = a 7 X 7 + Xs + X 9 + ^] 3 =10 a i X i , 

3) X = a 6 X 6 + arXr + X 8 +X 9 + J2]i 11 aiX i , 

4) X = a 5 X 5 + X 8 + X 9 + auXi + ai 3 X 13 , 

5) X = Y,l= 5 a ^+ X ^+ X a + a ^is, 

6) X = a 4 X 4 + X 9 + J2ll 10 OiXi, 

7) X = X 4 +a5X B + o 8 X 8 + X;" 10 OiX i , 

8) X — Xa + aeXtj + a$X$ + aioXio + arzXrz + 013X13, 

9) X — X4 + 07X7 + a$X 8 + aioXio + 012X12 + 013X13, 

10) X = X 4 + ELb a ^ + a 13Xl3, 

11) x = x 4 + a 7 x 7 + E" 10 "iXi, 

12) X = X 7 + ^ 10 aiX i , 

13) X = a3X 3 + J2ll 10 aiXi, 

U) X = X 4 + a 6 X 6 + J2Zu a > X i> 

15) X — X4 + 05X5 + 011X11 + 013X13, 

16) X = X 4 + Y,L 5 a * x * + 013X13, 
J 7; X = X 6 + X)" u OiXi, 

18) X = a 2 X 2 + ai Xi, 

19) X = a 2 X 2 + 012X12 + 013X13, 

20) X = X s + 011X11 + 013X13, 

21) X = a 2 X 2 + 03X3 + X 5 + 013X13, 
22; X = aiXi + Xn + 013X13, 

23) X = aiXi + 03X3 + X 6 + 013X13, 

24) X = aiXi + 02X2 + 03X3 + 013X13, 

25) X = Xg + aioXio + anXn + 013X13, 

26) X = a 3 X 3 + X10 + 011X11 + 013X13. 

Proof: F(sXi) : Q — > g defined byln Ad(exp(sXi)X) is a linear map, 
for i = 1, • • • , 13. The matrix of F(sXi), i = 1, ■ • • , 13, with respect 

to basis {Xi, • • • , X13} is 

1 0000000000" 

1 000 00 0000 

1 0000000000 

1000000000 

0100000000 

0010000000 

0001000000 

0000100000 

-s 0000010000 

s 0000001000 

000 000100 

O-sOOO 000010 

0000000001 
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and /13, respectively. Let A = cnXi, then 

F(s7.X 7 )oF(s 6 .X 6 )o-..oF(si.Xi) : X i-> 

(s3ai2 — S2ttio — siag + ai — S5S4CI8 + sr,a4 — S4as).Ai 

+ ( — S6S408 + S3«ll — S2CI9 + SiaiO +132+ Sea4 — S406).X2 

+(— S7S4O8 — S3i9 — s 2 au — s 1 a 12 + a3 + s 7 a 4 — s 4 a 7 ).X 3 (11) 
+ (a 4 - s 4 a s ).A 4 + (-s 6 a 10 - s 5 a g + s 5 a 8 + a 5 + s 7 ai 2 ).X 5 
+ ( — sea9 + S6a 8 + ssaio + a 6 + S7an).X 6 

+ (S708 — «6dll — Sodl2 + 0,7 — Sjat^.Xj 

+as-Xs + ag.Xg + aio-Aio + fln.Xn + ai 2 .Xr2 + 013. X13 

If as, and 09 7^ 0, and a 8 7^ 2a9, then we can make the coefficients of 
Xi, - ■ ■ ,Xj vanish. Scaling X if necessary, we can assume that as = 1. 
And X is reduced to Case 1. 

If as — ag / 0, and a w 7^ 0, then we can make the coefficients of 
Xi ■ ■ ■ ,Xe vanish. Scaling X if necessary, we can assume that os = 1. 
And X is reduced to Case 2. 

If as = ag / 0, aio = 0, and a\ 2 0, then we can make the coefficients 
of X\ , • • • , X5, and Aio vanish. Scaling X if necessary, we can assume that 
a8 = l- And X is reduced to Case 3. 

If as = ag ^ 0, aio = a\ 2 = 0, and an 0, then we can make the 
coefficients of Ai,---,A4, Xe, X7, Aio, and X12 vanish. Scaling X if 
necessary, we can assume that as = 1. And X is reduced to Case 4. 

If as = ag 7^ 0, aio = an = «12 = 0, then we can make the coefficients 
of Xi, ■ ■ ■ , X4, and X10, • ■ ■ , An vanish. Scaling X if necessary, we can 
assume that as = 1. And A is reduced to Case 5. 

If as =0, and ag 7^ 0, then we can make the coefficients of Xi, ■ ■ ■ , X3, 
and X5, ■ ■ ■ , Xs vanish. Scaling A if necessary, we can assume that ag — 1. 
And X is reduced to Case 6. 

If as 0, ag = 0, and an 7^ 0, then we can make the coefficients of 
Xi, • ■ ■ , A3, and Xq, Xi, and Xg vanish. And A is reduced to Case 7. 

If as 7^ 0, 09 = an = 0, and ai2 7^ 0, then we can make the coefficients 
of Ai, • • • , A3, and Xs, X7, Xg, and An vanish. And X is reduced to 
Case 8. 

If as 7^ 0, 09 = an = ai2 = 0, and aio 7^ 0, then we can make the 
coefficients of Xi, ■ ■ ■ , A3, and X$, X§, Xg, and An vanish. And X is 
reduced to Case 9. 
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If as / 0, and ag = aio = an = ai2 = 0, then we can make the 
coefficients of Ai, ■ ■ ■ , A3, and Xg, ■ ■ ■ , A12 vanish. And X is reduced to 
Case 10. 

If as = ag = 0, 04 =fc 0, and 010 7^ 0, then we can make the coefficients 
of X\, ■ ■ ■ , A3, A5, Xq, X$, and Xg vanish. Scaling X if necessary, we can 
assume that 04 = 1. And X is reduced to Case 11. 

If 04 = as = 09 = 0, aio 7^ 0, ai2 7^ 0, and a6ai2 + 07010 7^ asan, then 
we can make the coefficients of Ai, ■ ■ ■ ,X&, X$, and Xg vanish. Scaling 
X if necessary, we can assume that 07 = 1. And X is reduced to Case 12. 

If 04 = as = 09 = 0, aio 7^ 0, 012 7^ 0, and 00012 + 07010 = 05011, then 
we can make the coefficients of X\, X2, and X4, ■ ■ ■ ,Xg vanish. And X is 
reduced to Case 13. 

If as = ag = aio = 0, and 04 7^ 0, then we can make the coefficients of 
Xi, • • • , A3, A5, and Xf, ■ ■ ■ , Aio vanish. Scaling X if necessary, we can 
assume that 04 = 1. And X is reduced to Case 14. 

If as = ■ ■ ■ = aio = 0, 012 = 0, 04 / 0, and an / 0, then we can 
make the coefficients of X\, ■ ■ ■ , A3, and As, • • • , Aio vanish. Scaling A if 
necessary, we can assume that 04 = 1. And A is reduced to Case 15. 

If as = • ■ • = 012 = 0, and 04 7^ 0, then we can make the coefficients 
of Ai, • • • , A3, and As, • • • , Aio vanish. Scaling A if necessary, we can 
assume that 04 = 1. And A is reduced to Case 16. 

If 04 = as = 09 = aio = 0, 012 7^ 0, and a 5 an / 06012, then we can 
make the coefficients of Ai, • ■ ■ , A5, and A7, • • • , Aio vanish. Scaling A if 
necessary, we can assume that a6 = 1. And X is reduced to Case 17. 

If 04 = as = 09 = aio = 0, 012 7^ 0, asan = a6ai2, and 06 7^ 0, then 
we can make the coefficients of Ai, • • • , A5, and A7, • ■ ■ , Aio vanish. And 
X is reduced to Case 18. 

If 04 = oo = as = 09 = 010 = 0, 012 7^ 0, asan = 060i2, and ae 7^ 0, 
then we can make the coefficients of Ai, • ■ • , A5, and A7, • ■ ■ , Aio vanish. 
And X is reduced to Case 18. 

If 04 = a6 = as = • • • = on = 0, 012 7^ 0, 05011 = 060i2, 06 7^ 0, and 
a 5 / 0, then we can make the coefficients of Ai, and A3, • • ■ , An vanish. 
And X is reduced to Case 19. 

If 04 = OS = 09 = 010 = 012 = 0, 05 / 0, and an / 0, then we can 
make the coefficients of Ai, • • • , 04, Xe, ■ ■ ■ , Aio, and ai2 vanish. Scaling 
A if necessary, we can assume that 05 = 1. And A is reduced to Case 20. 

If 04 = as = • • • = 012 = 0, and 05 7^ 0, then we can make the 
coefficients of Ai, 04, and As, • • • , A12 vanish. Scaling A if necessary, we 
can assume that 05 = 1. And A is reduced to Case 21. 

If 04 = a 5 = as = 09 = aio = 012 = 0, and an / 0, then we can make 
the coefficients of A2, • • • , Aio and 012 vanish. Scaling A if necessary, we 
can assume that an = 1. And A is reduced to Case 22. 

If 04 = 05 = OS = • • • = 012 = 0, and ag 7^ 0, then we can make the 
coefficients of A2, 05, and A7, • • • , A12 vanish. Scaling A if necessary, we 
can assume that a& = 1. And A is reduced to Case 23. 

If 04 = 05 = oo = 0, and as = • • • = 012 = 0, then we can make the 
coefficients of A4, • ■ ■ , A12 vanish. And X is reduced to Case 24. 
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If an = as = ag = ai2 = 0, aio 7^ 0, and araio 7^ asan, then we 
can make the coefficients of X\, ■ ■ ■ , X4, Xe, ■ ■ ■ ,Xg, and a\2 — vanish. 
Scaling X if necessary, we can assume that as = 1. And X is reduced to 
Case 25. 

If C14 = as = ag = ai2 = 0, 07010 = asan, and aio 7^ 0, then we 
can make the coefficients of Xi, ■ ■ ■ , -X4, X@, ■ ■ ■ ,Xg, and 012 = vanish. 
Scaling X if necessary, we can assume that aio = 1- And X is reduced to 
Case 26. And the Theorem follows. □ 

According to our optimal system of one-dimensional subalgebras of the 
full symmetry algebra g, we need only find group-invariant solutions for 
26 one-parameter subgroups generated by X as Theorem 6. For example, 
as a direct consequence of Theorem 3 and case 1 of Theorem 6, we deduce 
the following conclusion: 

Conclusion 2 Let s, ci, • • • ,ce £ K with a, C2, and c$ =fc 0, and s be 

sufficiently small. Then, if 

u = U(t,x,y,z), v = V(t,x,y,z), w = W(t,x,y,z), 
p = P(t, x,y,z), q = Q(t, x, y, 2), 

is a solution of the system (1), so are the functions 

u — CiCjp ( cos(c3s) cos(c5s) — sin(c3s) sin(c4s) sin(css)) -U(t, x, y, z) 

— c i c 2~ S ( sin(c3s) cos(c 5 s) + cos(c3s) sin(c4s) cos(css)) -V(t, x, y, z) 

+CiC2" s cos(c4s) sm(css).W(t, x, y, z), 

v — CiCj 3 sin(c3«) cos(c4s).?7(t, x, y,z) 
+C1C2 s cos(c3«) cos(c4s).V(t, x, y, z) 
+c{c2 a . sin(c4s).W(t, x, y, z), 

w = —c\c^ s ( cos(c3s) sin(c 5 s) + sin(c3s) sin(c4s) cos(css)) -U(t, x, y, z) 
+ c i c 2 S ( sin(c3s) sin(c 5 s) — cos(c3s) sin(c4s) cos(css)) .V (t, x, y, z) 
+CiC2" S cos(c4s) cos(c5s).W(t, x, y, z), 

P = eg 8 ' .P(t,x,y,z), 

q = Ci c 2 c 6 .Q(t,x,y,z), 

where t = cf .t, and 

x = cl ( sin(css) sin(c4s) sin(c3«) + cos(css) cos(c3s)) .x 
+C2 ( sin(c3s) cos(css) — cos(c3s) sin(c4«) sin(c 5 s)) .y 
— c 3 2 cos(c4s) sin(c5s).2, 

y = — C2 sin(c3s) cos(c4s).a: + (C2S) cos(c3s) cos(c4s).y — C2 sin(c4s).2;, 

z = C2 ( cos(c 3 s) sin(c 5 s) — sin(c 3 s) sin(c 4 s) cos(c 5 s)).:r 
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+c 2 ( cos(c3s) sin(c4s) cos(css) + sin(c3s) sin(css)) .y 

+C2 C0S(C4S) COs(C5fi).2. 



The construction of the group- invariant solutions for each of the one- 
dimensional subgroups in the optimal system proceeds in the same fashion. 
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1. Introduction 

The equations governing the unsteady flow of a perfect polytropic gas in 
three space dimensions are 



where t is the time and x, y, z are the space coordinates; q(t,x,y,z) is the 
density, p(t,x,y,z) the pressure, u(t, x, y, z), v(t,x,y,z) and w(t, x,y,z) the 
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dw dw dw dw 1 dp 
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velocity components in the x, y and z directions, respectively, and 7 G R is the 
adiabatic index. These equations are called " three dimensional Euler equations 
of gas dynamics" [8] and section 6.3 of [9]. 

In this paper, we consider the equations governing the unsteady flow of a 
perfect polytropic gas in three space dimensions. The basic similarity reduc- 
tions for this system are performed. Reduced equations and exact solutions 
associated with the symmetries are obtained. We find an optimal system of 
one-dimensional sub-algebras for symmetry algebra of this system. 

In section 2 we find the full symmetry algebra of system (1). Chapter 3 is 
devoted to the group-invariant solutions to the system. The structure of full 
symmetry algebra illustrated in section. The last section deals with a optimal 
system of sub-algebras. 

2. Lie symmetry of the system 

The classical method for finding symmetry reductions of PDE is the Lie 
group method of infinitesimal transformations. To apply the classical method 
to (1), we consider the one-parameter Lie group of infinitesimal transformations 
in (t, x, y, z, u, v, p, q) given by 



t 


= t + s.^{t,x,y,z, 


u, 


v,p,q) + 0{s 2 ), 


X 


= x + s.&{t,x,y,z. 


u, 


V,p, q) + 0(s 2 ), 


y 


= y + s.£ 3 {t,x,y,z, 


u, 


V,p, q) + 0(s 2 ), 


z 


= z + s.£ 4 (t,x,y,z, 


u, 


v,p,q) + 0(s 2 ), 


u 


= u + s.rji(t,x,y, z. 


,11. 


,v,p,q) + 0(s 2 ), 


V 


= v + s.r] 2 {t 7 x,y,z : 


u. 


v,p,q) + 0(s 2 ), 


Ul 


= w + s.r] 3 (t 7 x,y,z 


,u 


,v,p,q) + 0(s 2 ), 


p 


= p + s.rj 4 (t,x,y,z, 


u, 


v,p,q) + 0(s 2 ), 


Q 


= Q + s-V5(t,x,y,z, 


u, 


V,p, q) + 0(s 2 ). 



where s is the group parameter. One requires that this transformation leaves 
the sub-manifold S invariant. Let 

X = £1 d t + £2 d x + £3 d y + £4 d z + ni d u + 772 d v + t? 3 d w + r\\ d p + 775 d q 

be the corresponding infinitesimal transformation; where coefficients are real- 
valued C 00 — functions of t, x, y, z, u, v, w,p, q. The infinitesimal criterion of in- 
variance (1) will thus involve t, x, y, z, u, v, w, p, q and the derivatives of u, v, w,p, q 
with respect to t, x, y, z, as well as £j, i = 1, 2, 3, 4, and rjj, j — 1, 2, 3, 4, 5, and 
their partial derivatives. After eliminating any dependencies among the deriva- 
tives of the m's caused by the system itself, we can then equate the coefficients 
of the remaining unconstrained partial derivatives of u to zero. This will result 
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in a large number of elementary partial differential equations for the coefficient 
functions r]j of the infinitesimal generator, called the determining equations 
for the symmetry group of the given system. In this case, we find the determin- 
ing equations for the symmetry group of the system (1) to be the following: 



Cm 2 


= £i,x 


- €l,V 


= Sl,z 


= £l,u 








= £i, P = 0, 




— £2,u 


— &,v 


— &,w 


= 6,9 


= 6,p 


= 6,t 2 


= &,ty 


= &,y 2 = 0, 




= £3,11 


— &,v 


— &,w 


= 6,9 




= 6,t 2 


= o, 




U,z 




= U,q 


= ^4,p 


= 0, 











£,i,t 2 — £4,te — £4,tj/ — £4, z 2 — Z.i,xy — ^4,y 2 — 0, 

£,3,x = £2,1/1 &,z + U,y = 0, = £4,1*1 £2,2 + £,4,,x = 0, 

and 

Vl = _u £l,t + 6,t + U&.i; - w£ 4j:E , 
^72 = -w£,4,y + &,t - U&,y ~ u£l,t, 

m = u U,x + U,t + v£ 4 , v - ty£i,t, (4) 

V5,t = V5,x = V5,y = V5,z = V5,u = V5,v = V5,w = V5,q = 0, 
V5 = PVb.p- 

First, equations (3) require that £jS are just 
£1 = a-it + o 2 , 

£2 = ^3* + 04a; + a 5 y + a & z + a 7 , 

£3 = a$t - a 5 x + a A y + a 9 z + a w , (5) 
£4 = ant — ciqx — agy + a^z + 012, 

where a\, ■ ■ ■ , a 12 are arbitrary constants. Then, equations (4) and (5) requires 
that 

f]i = (04 - ai)u + a 5 v - a 6 w + 03, 
f]2 = -ar-,u + (a 4 - a\)v + a^w + a 8 , 

V3 = -a&u - a 9 v + (a 4 - ai)w + an, (6) 
774 = (2ai - 2a 4 + a 13 )q, 
V5 = a 13 p, 

where 013 is an arbitrary constant. Therefore 
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Theorem 1.. The Lie algebra of infinitesimal symmetries of the system (1) is 
spanned by the 13 vector fields 



x 1 


= d x , 




x 2 


= 9y, 


(shifts of origin) 


x 3 


= d z , 




x 4 


= d t , 




x 5 


= td x + d u , 




x 6 


= tdy + 8y, 


( uniformly frame motions ) 


x 7 


= td z + d w , 




X & 


= t dt — u d u — v d v — w d w 




x 9 


= x d x + y d y + z d z + u d u 


+ v d v + w d w — 2r d r , 



(7) 



X w = y d x - x d y + v d u - u d v , 

Xn = —zd y + y d z — w d v + vd w , (rotations of reference frame) 
X 12 = -z d x + x d z - w d u + u d w , 



Xi 3 ^qd q +pd p , 

These infinitesimal symmetries will generate a Lie algebra g over the field of 
real or complex numbers. The commutator table of Lie algebra g for (1) is given 
below, where the entry in the ith row and jth column is defined as [Xi,Xj] = 
Xi.Xj - Xj.Xi, i,j = !,■■■, 13. 





Xi 


x 2 


x 3 


X 4 


x 5 


x 6 


x 7 


Xg Xg 


X10 


Xn 


Xi 2 


Xl 3 

























Xj 


-X a 





x 3 





x 2 























x 2 


Xi 


x 3 








x 3 























x 3 





-x 2 


-Xi 





Xa 














Xi 


x 2 


x 3 


Xa 














x 5 











-Xi 











— 


-x 6 





X 7 





x 6 











-x 2 











— Xq Xq 


x 5 


X r 








x 7 











-x 3 











— X7 X7 





-X 6 


-x B 





x 8 











-Xa 


x 5 


x 6 


X 7 

















x 9 




-x 2 


-x 3 





-x 5 


-x 6 


-X 7 

















Xu) 


x 2 










x 6 


-X t 











— Xi 2 


X n 





Xn 





-x 3 


x 2 








-X T 


x 6 





Xi 2 





— X10 





Xl2 


-X 3 





Xi 





-X r 





x 5 





-Xu 


X10 








Xl 3 







































3. Reduction of the system 

Now, we look at the group-invariant solutions to the system (1). 
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The system (1) is a sub-manifold M of jet space J 2 (R 4 ,R 5 ) defined by the 
following equations 

u t + uu x + vuy + wu z +p x /q = 0, 
v t + uv x + vv y + wv z + py/q = 0, 

w t + uw x + vwy + ww z +p z /q = 0, (8) 
qt + q (u x +v y + w z ) + uq x + vq y + wq z = 0, 

Pt + IP {U X +Vy+ W Z ) + UP X + VPy + WP Z = 0. 

Doing as section 3.1 of [6] and find (in a sense) the most general group- 
invariant solutions to the system (1). 

Theorem 2.. The one-parameter groups gi(s) : M — > M generated by the Xi, 
i = 1, • • • , 13 are given in the following table: 

91 ■ (t,x + s,y,z,u,v,w,p,q) 

92 ■ (t,x,y + s,z,u,v,w,p,q) 

93 ■ (t,x,y,z + s,u,v,w,p,q) 

94 ■ (t + s,x,y,z,u,v,w,p,q) 

95 ■ (t,x + st,y,z,u + s,v,w,p,q) 
9e ■ (t,x,y + st,z,u,v + s,w,p,q) 

97 ■ (t,x,y,z + st,u,v,w + s,p,q) (9) 
,g 8 : (e s i, x, y, z, e~ s u, e' s v, e~ s w,p, e 2s q) 
gg : (t,e s x,e s y,e s z,e s u,e s v,e s w,p,e~ 2s q) 
gio : (t, xcos s + y sins, y cos s — a; sins, z, 

v sin s + u cos s, v cos s — u sin s, w, p, q) 
gn : (t, x, y cos s — z sin s, z cos s + y sin s, 

u, v cos s — w sin s, w cos s + v sin s, p, q) 
gu : (i, xcos s — zsins, y, resins + zcoss, 

u cos s — w sin s, v, w cos s + u sin s, p, q) 
313 : (t,x,y,z,u,v,w,e s p,e s q) 

Where, entries give the transformed point 

cxp(sX l )(t,x,y,z,u,v,w,p,q) = (t,x,y,z,u,v,w,p,q). 

Since each group gi, i = 1, ■ • • , 13, is a symmetry group of (1), then 
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Theorem 3.. Let i = 1, • • • , 13, s G R, and A G K - {0}. 7/ 

u = U(t,x,y,z), v = V(t,x,y,z), w = W(t,x,y,z), 
p = P(t, x, y, z), q = Q(t, x, y, z), 

is a solution of the system (1), so are the functions 

=U®(t,x,y,z), v® = V®(t,x,y,z), w w = W®(t,x,y,z), 
P W=pW{t,x,y,z), q®=Q®{t,x,y,z), 

where 





= U(t,x + s,y,z), 




= ?/(t,a:,2/ + s,z), 


V W 


= V(t,x + s,y,z), 


„(2) 


= V"(t,a;,j/ + s,z), 


W W 


= W(t,x + s,y,z), 


U,(2) 


= W(t,x,y + s,z) 


p(i) 


= P(t,x + s,y,z), 


p (2) 


= P(t,x,y + s,z), 


q W 


= Q(t,x + s,y,z). 


g (2) 


= Q(t,x,y + s,z). 



M ( 3 ) = U(t,x,y, z + s), 
v^l = V(t,x,y,z + s), 
w ( 3 ) = W(t,x,y,z + s), 
p( 3 ) = P(t,x,y,z + s), 
g (3) = Q(t,x,y,z + s). 





= f/(t + s 


,x,y- 


■*) 


V W 


= V(t + s. 


• x,y., 


z). 


W W 


= W(t + . 


s,x,y,z 


p (4) 


= + 


x,y, 


z), 


g (4) 


= Q(t + s, 


x,y, 


z). 



w ( 5 ) = {/(f, a; + st,y, z) - s, 

= V(t,x + st,y,z), 
u,(5) = W(t,x + st,y,z), 
p< 5 ) = P(t,x + st,y,z), 
q {5) = Q(t,x + st,y,z). 



u^l = U(t,x,y + st,z), 
w ( 6 ) = V(t,x,y + st,z) -s, 
w ( 6 ) = W(t,x,y + st,z), 
p(«) = P(t,a;,y + «*,«), 
g (6) = Q(t,x,y + st,z). 



u ( 7 ) = U(t,x,y,z + st), 
V <J) = V(t,x,y,z + st), 

= W(t,x,y,z + st) - s, 
P W = P(t,x,y,z + st), 
<? (7) = Q(t,x,y,z + st). 



= X.U(X.t,x,y,z), 
«(8) = A.V(A.t,a;,y,«), 

= X.W(X.t,x,y,z), 
p^ = P(X.t,x,y,z), 
q (8) = ^.Q(X.t,x,y,z). 



u ( V = \.U(t,X.x,X.y,X.z), 
w (9) = i.V(i,A.a;,A.y,A.z), 
w (9) = i.V^(t,A.x, A.y, A.z), 

= P{t,X.x,X.y,X.z), 
q^ = X 2 .Q(t,X.x,X.y,X.z). 



U ( 13 ) = U(t,x,y,z), 

p( 13 ) - A.P(t, 

g ( 13 ) = \.Q(t,x,y,z). 
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u' 10 -* = cos s.U(t, xcoss + y sins, ycoss — xsins, z) 

— sin s.V(t, x cos s + y sin s, y cos s — x sin s, z) , 
t/ 10 ' = sins. U(t, xcoss + y sins, ycoss — xsins, z) 

+ cos s.V^i, xcoss + y sins, ycoss — xsins, z), 
u/ 10 ) — IT (i, xcoss + ysins, ycoss — xsins, z), 
p(io) _ p(t, xcoss + ysins, ycoss — xsins, z), 
qt 10 ' 1 = xcoss + y sins, ycoss — xsins, z). 

u( n ) = U (t,x,y cos s — zsins, zcoss + y sins), 
i/ 11 ) = cos s.y(t,x, ycoss — zsins, zcoss + ysins) 

+ sins. IT (i, x, y cos s — z sins, z cos s + y sin s) , 
u/ 11 ) = cos s.W(t, x, y cos s — z sin s, z cos s + y sin s) 

— sins.V r (t, x, ycoss — zsins, zcoss + ysins), 
p( n ) _ p(t,x, ycoss — zsins, zcoss + ysins), 

</ n ) = Q(t,x, ycoss — zsins, zcoss + ysins). 

w' 12 ) = cos s.U(t, xcoss — zsins, y, xsins + zcoss) 

+ sin s. IT (i, xcoss — z sins, y, xsins + z coss), 
v ( 12 ) _ V(t, xcoss — zsins, y, xsins + z coss), 
u/ 12 ) = coss. W(t, xcoss — zsins, y, xsins + z coss) 

— sin s.U(t, xcoss — z sin s, y, x sins + zcoss), 
p( 12 ) = P(t, xcoss — zsins, y, xsins + z coss), 

</ 12 ^ = Q(t, xcoss — zsins, y, xsins + zcoss). 



4. Structure of Lie algebra g 

In this section, we determine the structure of full symmetry algebra g. of 
system (1). 

The center 3 of g is Span R {Xi 3 }. Therefore, the quotient algebra gi := g/3 
is Span R {Yi, • • • , Y12}; where Yi := Xi +3, i = 1, ■ • • , 12. The commutator table 
of Lie algebra gi is given below, where the entry in the ith row and jth column 
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is defined as [Y^Yj] = Y l .Y ] - Yj.Yi, i,j = !,■■■, 12. 





Yi 


Y 2 


Y 3 


Y A 


Y 5 


Y 6 


Y 7 


Yr 


Y 9 


Yio 


Yu 


Yi2 




























Yi 


-Y 2 





Y 3 


Y 2 


























Y 2 


Yi 


Y 3 





Y 3 


























Y 3 





-Y 2 


-Yi 


Y 4 














Yi 


Y 2 


Y 3 


Yi 














Y s 











-Yi 











-Y B Y 5 


-Y 6 





Y 7 


Y 6 











-Y 2 











-Y 6 Y 6 


Ys 


Y 7 





Y 7 











-Y 3 











-Y r 


Y 7 





-Ye 


-Yb 


Y 8 











-Yt 


Y 5 


Y 6 


Y 7 

















Y 9 


-Yi 


-Y 2 


-Y 3 





-Y 5 


-Y 6 


-Y, 

















Y w 


Y 2 


-Yi 








Ye 


-Y B 














-Yu 


Yu 


Yu 





-Y 3 


Y 2 








-Yi 


Y 6 








Y 12 





-Y w 


Yl 2 


-Y 3 





Yi 





-Y r 





Y 5 








-Yu 


Y w 






The algebra gi is not solvable, because 

si^ ■= [fli.fli] =Span M {Y 1 ,---,Y 10 }, and flf ) := [fl^, £1^] = fl^. 
Lie algebra gi admits a Levi-decomposition into the semi-direct sum 

01 = t K 5, 

where r = Span R {Y l7 • • • , Y 9 } is the radical of 0i (the largest solvable ideal 
contained in jji) and s = Span R {Yi , Yn, Yi 2 } is a semi-simple subalgebra 
of fli. The radical r is uniquely defined but the semi-simple subalgebra s is 
not. s as an algebra is isomorphic to 3— dimensional special orthogonal algebra 
so(3) := {A e Mat(3 x 3;R) | A* = -A}, which is a simple algebra. 
The radical r is solvable, with and 

t« =Span R {Yi,...,Y 7 }, t< 2 > = Span R {Yi,Y 2 ,Y 3 }, t< 3 > = {0}. 

It is semi-direct sum r = r^ 1 ) K 02, where 02 = Span R {Ys, Yg} is isomorphic to 
the two dimensional Abelian Lie algebra K 2 ; i.e. 

r = r (1) k R 2 . 

r' 1 ) is semi-direct sum t^ 1 ) = r* 2 ^ K 3 , where 3 = Span R {Y 4 , • • • , Y 7 } is 
isomorphic to the four dimensional Abelian Lie algebra R 4 , and is isomorphic 
to the three dimensional Abelian Lie algebra R 3 ; i.e. 

r (1) = R 3 k R 4 . 

Theorem 4-- The symmetry algebra of system (1) have the following structure 
= M x (((R 3 x R 4 ) k M 2 ) x so(3)). 
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Theorem 5.. There is a normal Lie-subalgebra of symmetry Lie-group of system 
(1) which is isomorphic to the Galilean group o/R 4 : 

r/i o ,\ 

Gal(4,R)=M v R y R G 0(3, R), s e R, and y, v e R 

{ \ o o i y 

Therefore, the system (1) is invariant up to Galilean motions of space-times R 4 
(see [3]). 

Proof. Let b = Span R {Xi, ■ • • , X7, X10, Xn, X12}. Lie algebra structure of b 
is as flo[(3), the lie algebra of Galilean group Gal(4,R), and b an ideal of g. 
Therefore, there is a Lie-subgroup of G such that its Lie-algebra is b, by the 
Theorem 2.53 of Olver [6]. □ 

Conclusion 1.. q is semi-direct sum of 

fl0l(4) = Span R {X 1 , • • • , X 7 , X 10 , X n , X 12 } 

and 3— dimensional Abelian algebra R 3 = Span{X 8 , Xg, X13}. 

5. Optimal system of sub-algebras 

As is well known, the Lie group theoretic method plays an important role 
in finding exact solutions and performing symmetry reductions of differential 
equations. Since any linear combination of infinitesimal generators is also an 
infinitesimal generator, there are always infinitely many different symmetry sub- 
groups for the differential equation. So, a mean of determining which subgroups 
would give essentially different types of solutions is necessary and significant for 
a complete understanding of the invariant solutions. As any transformation in 
the full symmetry group maps a solution to another solution, it is sufficient to 
find invariant solutions which are not related by transformations in the full sym- 
metry group, this has led to the concept of an optimal system [7] . The problem 
of finding an optimal system of subgroups is equivalent to that of finding an 
optimal system of subalgebras. For one-dimensional subalgebras, this classifica- 
tion problem is essentially the same as the problem of classifying the orbits of 
the adjoint representation. This problem is attacked by the naive approach of 
taking a general element in the Lie algebra and subjecting it to various adjoint 
transformations so as to simplify it as much as possible. The idea of using the 
adjoint representation to classify group-invariant solutions was due to [7] and 
[5]. 

The adjoint action is given by the Lie series 

Ad(exp(sX i )X j ) = Xj - s[Xi,Xj] + ^-[X i; [X h Xj]] - • ••, (10) 

where [Xi,Xj] is the commutator for the Lie algebra, s is a parameter, and 
i, j = 1, • • • , 13. We can write the adjoint action for the Lie algebra g, and show 
that 
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Theorem 6.. A one- dimensional optimal system of (1) is given 



1 ) X 

1J A 




To -1- V 13 n X - 


Z) A 




07A7 + A 8 + Ag + 2-,i=w a « A i; 


oj A 




a 6 A 6 + a 7 A.7 + A 8 + Ag + 2^i=ll a i A i; 


), ) X 




05^.5 ^ ^8 ^9 ^ "-ll^z 043W13, 


5 J X 




V 7 * a,-X + Xs + Xa + ffli qXi •} 


6) X 




^ /i a — 1— — 1— ^ ^ . _ „ fl a \: 
u '4 yi -4 \ -'•-y / y ^ — 10 * * 7 






Y\ _|_ -u no Xa -4- V 13 n ■ 

i "5^5 * "8^8 * z_^i— 10 u t^-tf 


8) X 




-X*4 H - 0-6-^6 H - GgXg ~h Qio-^lO H~ ^12^12 H~ &13-^13? 


9) X 




4 I ^11 > o O 1 1U1U 1 LZ 1 Z > 1 LO 1 


10) X 




Xa + K CLjXj ~h <2l 3-^13 5 


11) X 




Xa -\- CL7X7 H - ^ ■ -1 CLiXi , 

1 1 I I 1 / jq — \\j t i ? 


19) X 




y_ 1 v^ 13 „.v-. 
A7 -|- z^i=io 8 ' A 'i 


10 J A 






U) x 






15) X 




X4 + Asa's + flulu + 013X13, 


16) X 


— 


X4 + Z)I=5 a *-^» + a 13-Xl3> 


17) X 


= 


^6 + 5^1=11 a i^i; 


18) X 




02X2 + Si=ll Oi^i; 


19) X 




^2X2 + 012X12 + 013X13, 


20) X 




X5 + aiiXn + 013X13, 


21) X 




a 2 X 2 + a 3 X 3 + X 5 + 013X13, 


22) X 




aiXi + Xn + 013X13, 


23) X 




aiXi + a 3 X 3 + X 6 + 013X13, 


24) x 




aiXi + a 2 X 2 + a 3 X 3 + 013X13, 


25) X 




X5 + oiqXio + OiiXn + 013X13, 
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26) X = a 3 X 3 + X 10 + auXu + a 13 X 13 . 



Proof: F(sXi) : g — > g defined by X h Ad(cxp(sX;)X) is a linear map, for 
i = 1, • ■ • , 13. The matrix M s x t of F(sXi), i = 1, • ■ ■ , 13, with respect to basis 
{X 1 ,---,X 13 }is 
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and -Z13, respectively. Let X = ^,=i a iXi, then 

%.X 7 )o%.X 6 )o-.-oF(s 1 .X 1 ) : X i ► 

(s 3 a 12 - s 2 aio - «lfl9 + Ol - S5 s 4fl8 + s 5 a 4 - s 4 a 5 ).X! 

+ (-s e s i a 8 + s^an - s 2 a g + sia w + a 2 + s e a 4 - s i a 6 ).X 2 

+ (-s 7 s i a s - s 3 a 9 - s 2 an - S1O12 + a3 + s 7 a 4 - s 4 a 7 ).X 3 (11) 

+ (a 4 - s 4 a 8 ).X 4 + (-S6O10 - 5509 + s 5 a s + a 5 + s 7 ai 2 ).X 5 

+ (-s 6 a 9 + s 6 a 8 + s 5 aio + a 6 + S7an).X 6 

+ (s7a 8 - s & an - s 5 a 12 + a 7 - s 7 a 9 ).X 7 

+a 8 .A 8 + ag.Xg + aiQ.X\Q + ffln-Xii + a\ 2 .X\ 2 + 013. X13 

If as, and ag 7^ 0, and a 8 7^ 2ag, then we can make the coefficients of 
X\, ■ ■ ■ , X 7 vanish. Scaling X if necessary, we can assume that as = 1. And X 
is reduced to Case 1. 

If a 8 = ag =/= 0, and aio 7^ 0, then we can make the coefficients of X\ • • • , Xq 
vanish. Scaling X if necessary, we can assume that a 8 = 1. And X is reduced 
to Case 2. 

If a 8 = a 9 7^ 0, a 10 = 0, and a 12 7^ 0, then we can make the coefficients of 
Xi, ■ ■ ■ , X 5 , and X w vanish. Scaling X if necessary, we can assume that a 8 = 1. 
And X is reduced to Case 3. 

If a 8 = ag 7^ 0, aio = a i2 = 0, and an 7^ 0, then we can make the coefficients 
of X\,- ■ ■ ,X±, Xq, X 7 , X\o, and X 12 vanish. Scaling X if necessary, we can 
assume that a 8 = 1. And X is reduced to Case 4. 

If as = «9 7^ 0, aio = an = ai2 = 0, then we can make the coefficients of 
Xi, ■ ■ ■ , X4, and X\o, • • • , Xn vanish. Scaling X if necessary, we can assume 
that as = 1. And X is reduced to Case 5. 

If as — 0, and a 9 7^ 0, then we can make the coefficients of X\, ■ ■ ■ , X3, and 
X 5 , ■ ■ ■ , X s vanish. Scaling X if necessary, we can assume that a 9 — 1. And X 
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is reduced to Case 6. 

If as 7^ 0, ag = 0, and an ^ 0, then we can make the coefficients of 
X\, ■ ■ ■ , A3, and Xq, A7, and Ag vanish. And X is reduced to Case 7. 

If as 7^ 0, ag = an = 0, and a\2 7^ 0, then we can make the coefficients of 
Xi, ■ ■ ■ , Xs, and A 5 , X-j, Ag, and An vanish. And X is reduced to Case 8. 

If ag 7^ 0, ag = an = a i2 = 0, and aio =/= 0, then we can make the coefficients 
of Xi, ■ ■ ■ , A 3 , and X 5 , A 6 , Ag, and An vanish. And X is reduced to Case 9. 

If ag ^ 0, and a 9 = a w = an = a i2 = 0, then we can make the coefficients 
of X\, ■ ■ ■ , A3, and Ag, • • • , A12 vanish. And A is reduced to Case 10. 

If as = ag = 0, 04 7^ 0, and aio 7^ 0, then we can make the coefficients of 
Ai, • • • , A 3 , A 5 , A 6 , Ag, and Ag vanish. Scaling A if necessary, we can assume 
that 04 = 1. And A is reduced to Case 11. 

If a4 = as = ag = 0, aio 7^ 0, ai2 0, and a6ai2 + 07010 7^ a^an, then 
we can make the coefficients of Ai, • • • , Xq, Xs, and Ag vanish. Scaling A if 
necessary, we can assume that 07 = 1. And A is reduced to Case 12. 

If 04 = as = ag = 0, aio 7^ 0, ai 2 7^ 0, and 06012 + 07010 = asan, then we 
can make the coefficients of Ai, A 2 , and A4, • • • , Ag vanish. And A is reduced 
to Case 13. 

If as = ag = aio = 0, and 04 7^ 0, then we can make the coefficients of 
Ai, • • • , A 3 , A 5 , and X?, ■ ■ ■ , A i0 vanish. Scaling A if necessary, we can assume 
that 04 = 1. And A is reduced to Case 14. 

If as = • • • = aio = 0, 012 = 0, 04 7^ 0, and an 7^ 0, then we can make the 
coefficients of Ai, • • • , A 3 , and A 8 , • • • , A i0 vanish. Scaling A if necessary, we 
can assume that 04 = 1. And A is reduced to Case 15. 

If as — ■ ■ ■ = oi 2 = 0, and 04 7^ 0, then we can make the coefficients of 
Ai,---,A 3 , and A 8 ,---,Ai vanish. Scaling A if necessary, we can assume 
that 04 = 1. And A is reduced to Case 16. 

If 04 = as = ag = aio = 0, ai 2 7^ 0, and 05011 7^ agai2, then we can make 
the coefficients of Ai, • • • , A 5 , and A7, • • • , Ai vanish. Scaling A if necessary, 
we can assume that a§ = 1. And A is reduced to Case 17. 

If 04 = as = ag = aio = 0, 012 7*= 0, a^an = a$ai2, and ae 7^ 0, then we can 
make the coefficients of Ai, • • • , A 5 , and A7, • • • , Aio vanish. And A is reduced 
to Case 18. 

If 04 = ag = a§ = ag = aio = 0, ai 2 7^ 0, asan = agai2, and a^ 7^ 0, then 
we can make the coefficients of Ai, • • • , A5, and A7, • • • , Aio vanish. And A is 
reduced to Case 18. 

If 04 = ae = as = ■ ■ ■ = an = 0, ai 2 7^ 0, a 5 an = a 6 ai 2 , a 6 7^ 0, and a 5 7^ 0, 
then we can make the coefficients of X\, and A3, • • • , An vanish. And A is 
reduced to Case 19. 

If 04 = ag = ag = aio = 012 = 0, 05 7^ 0, and an 7^ 0, then we can make the 
coefficients of Ai, • • • , 04, Xq, • • • , Aio, an d ai2 vanish. Scaling A if necessary, 
we can assume that a 5 = 1. And A is reduced to Case 20. 

If a 4 = ag = • • • = a i2 = 0, and a 5 7^ 0, then we can make the coefficients 
of Ai, a 4 , and Ag, • • • , A i2 vanish. Scaling A if necessary, we can assume that 
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as = 1. And X is reduced to Case 21. 

If 04 = 05 = as = ag = aio = ai2 = 0, and an 7^ 0, then we can make 
the coefficients of X2, ■ ■ ■ ,X\o and 012 vanish. Scaling X if necessary, we can 
assume that on = 1. And X is reduced to Case 22. 

If 0,4 = 05 = ag = ■ ■ ■ = ai2 = 0, and a 6 7^ 0, then we can make the 
coefficients of X 2 , 05, and X?, ■ ■ ■ ,X 12 vanish. Scaling X if necessary, we can 
assume that a 6 ~ 1. And X is reduced to Case 23. 

If a 4 = a 5 = a 6 = 0, and a 8 = • • • = a 12 = 0, then we can make the 
coefficients of X4, • • • , X12 vanish. And X is reduced to Case 24. 

If 04 = as = ag = ai2 = 0, aio 7^ 0, and 07010 7^ 05011, then we can make 
the coefficients of X\, ■ ■ ■ , X4, Xq, • • • , Xg, and 012 = vanish. Scaling X if 
necessary, we can assume that 05 = 1. And X is reduced to Case 25. 

If 04 = as = ag = ai2 = 0, 07010 = 05011, and oio 7^ 0, then we can make 
the coefficients of X\, ■ ■ ■ , X4, X$, • • • , Xg, and 012 = vanish. Scaling X if 
necessary, we can assume that Oio = 1. And X is reduced to Case 26. And the 
Theorem follows. □ 

According to our optimal system of one-dimensional subalgebras of the full 
symmetry algebra g, we need only find group- invariant solutions for 26 one- 
parameter subgroups generated by X as Theorem 6. For example, as a direct 
consequence of Theorem 3 and case 1 of Theorem 6, we deduce the following 
conclusion: 

Conclusion 2. Let s, Ci, • ■ • ,c% el with c\, C2, and cq ^ 0, and s be sufficiently 
small. Then, if 

u = U(t,x,y,z), v = V(t,x,y,z), w = W(t,x,y,z), 
P = P{t, x,y,z), q = Q(t, x, y, z), 

is a solution of the system (1), so are the functions 

u = clc2~ B ( cos(c3s) cos(c5s) — sin(c3s) sin(c4s) sin(c5s)) -U(t, x,y, z) 
— c\c2~ s ( sin(c3s) cos(c5s) + cos(c3s) sin(c4s) cos(c5s)) -V(t, x, y, z) 
+c{ c 2 s cos(c4s) sm(c5s).W(t, x, y, z), 

v = c\c2 S sin(c3s) cos(c4s).U(t, x, y, z) 
+c\c 2 s cos(c3s) cos(c4s).V(t, x, y, z) 
+c\c2~ s . sm(c4s).W(t, x, y, z), 

w = —c\c2 S ( cos(c3s) sin(c5s) + sin(c3s) sin(c4s) cos(css)) -U(t, x, y 7 z) 
+ c i c 2^ S ( sin(c3s) sin(c5s) — cos(c3s) sin(c4s) cos(c5s)) -V(t, x, y, z) 

+C S l C2 S C0S(C4S) COs(C5S).W(f, X, X), z), 

P = CQ S .P(i,x,y,z), 
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q = c 1 2s cl s c 6 s .Q(t,x,y,z), 

where t = c\.t, and 

x = c s 2 ( sin(c5s) sin(c4s) sin(c 3 s) + cos(css) cos(css)).x 
+C2( sin(c3s) cos(c5s) — cos(c3s) sin(c4s) sin(c 5 s)).j/ 
— c\ cos(c4s) sin(c5s).z, 

y = — sin(c3s) cos(c4s).a; + (c2s) cos(c3s) cos(c4s).y — C2 sin(c4s).z, 

z — c s 2 ( cos(c 3 s) sin(c 5 s) — sin(c 3 s) sin(c4s) cos(c5s)) .a; 

+C2 ( cos(c 3 s) sin(c4s) cos(c5s) + sin(c 3 s) sin(css)) .y 

+C 2 C0S(C4S) C0S(C5S).Z. 

The construction of the group-invariant solutions for each of the one-dimensional 
subgroups in the optimal system proceeds in the same fashion. 
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Abstract 

In this paper, the equations governing the unsteady flow of a perfect 
polytropic gas in three space dimensions are considered. The basic sim- 
ilarity reductions for this system are performed. Reduced equations and 
exact solutions associated with the symmetries are obtained. This results 
is used to give the invariance of system up to Galilean motions of space- 
times K 4 . Then, an optimal system of one-dimensional sub-algebras for 
symmetry algebra of this system is given. 

Keywords: Euler equations, Lie group of transformations, Symmetry 
generators, Similarity solution, optimal system of Lie sub-algebras. 
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1 Introduction 

The equations governing the unsteady flow of a perfect polytropic gas in 
three space dimensions are 

du du du du 1 dp 

dt dx dy dz q dx 

dv dv ^dv ^dv 1 dp _ ^ 

dt dx dy dz q dy 

dw dw dw dw 1 dp , . 

dt dx dy dz q dz ' 

dq ( du dv dw \ dq dq ^ dq _ ^ 

dt ^\dx dy dz ) dx dy dz 

dp / du dv dw \ dp ^ ^ dp ^ dp _ ^ 

dt ^ v dx dy dz ) dx dy dz 

where t is the time and x, y, z are the space coordinates; q(t, x, y, z) is the 
density, p(t,x,y,z) the pressure, u(t,x,y,z), v(t,x,y,z) and w(t,x,y, z) 
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the velocity components in the x, y and z directions, respectively, and 7 £ 
M is the adiabatic index. These equations are called "three dimensional 
Euler equations of gas dynamics" [8] and section 6.3 of [9] . 

In this paper, we consider the equations governing the unsteady flow 
of a perfect polytropic gas in three space dimensions. The basic similarity 
reductions for this system are performed. Reduced equations and exact 
solutions associated with the symmetries are obtained. We find an opti- 
mal system of one-dimensional sub-algebras for symmetry algebra of this 
system. 

In section 2 we find the full symmetry algebra of system (1). Chapter 
3 is devoted to the group-invariant solutions to the system. The structure 
of full symmetry algebra illustrated in section. The last section deals with 
a optimal system of sub-algebras. 

2 Lie symmetry of the system 

The classical method for finding symmetry reductions of PDE is the Lie 
group method of infinitesimal transformations. To apply the classical 
method to (I), we consider the one-parameter Lie group of infinitesimal 
transformations in (t,x,y,z,u,v,p,q) given by 



t 


= t + s.&(t,x,y,z, 


u, 


v,p,q) + 0(s 2 ), 


X 


= x + s.&(t,x,y,z, 


u, 


v,p,q) + 0(s 2 ), 


V 


= y + s.£ 3 (t,x,y,z, 


u, 


v,p,q) + 0(s 2 ), 


z 


= z + s.£ 4 (t,x,y,z, 


u, 


v,p,q) +0(s 2 ), 


u 


— u + s.r/i (t, x, y, z. 


, u. 


.v,p,q) + 0(s 2 ), 


V 


= v + s.rj 2 (t,x,y,z. 


u, 


v,p,q) + 0(s 2 ), 


w 


= w + s.r/ s (t,x,y,z 


,u 


,v,p,q) + 0(s 2 ), 


p 


= p + s.f] 4 (t,x,y,z, 


u, 


v,p,q) + 0(s 2 ), 


q 


= q + s.r) 5 (t,x,y,z, 


u, 


v,p,q) + 0(s 2 ). 



where s is the group parameter. One requires that this transformation 
leaves the sub-manifold S invariant. Let 

X = £1 d t + £2 d x + £3 d y + £4 d z + ??i d u + ??2 d v + m 9 W +V4.d p + r) 5 d q 

be the corresponding infinitesimal transformation; where coefficients are 
real-valued C°° —functions of t, x, y, z,u,v, w, p, q. The infinitesimal crite- 
rion of invariance (1) will thus involve t,x,y, z, u, v, w, p, q and the deriva- 
tives of u,v,w,p,q with respect to t,x,y,z, as well as i = 1,2,3,4, and 
Tjj, j = 1, 2, 3, 4, 5, and their partial derivatives. After eliminating any de- 
pendencies among the derivatives of the m's caused by the system itself, 
we can then equate the coefficients of the remaining unconstrained partial 
derivatives of u to zero. This will result in a large number of elementary 
partial differential equations for the coefficient functions £», r/j of the in- 
finitesimal generator, called the determining equations for the symmetry 
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group of the given system. In this case, we find the determining equations 
for the symmetry group of the system (1) to be the following: 

£i,t 2 = £i,x = £1,1/ = £i,z = = £i,« = £1,10 = £1,9 = £i, P = 0, 

£2,2; = £,2,u = &,v = £2,™ = £2,9 = £,2,p = £ 2 ,t 2 = £2,ty = £2, y 2 = 0, 

^3, y = £3,11 = £3,1) = £3,™ = £3,9 = £3,p = C3,t 2 = 0, (3) 

£4,2 = £,A,w = £,A,q = £4,p = 0, 

?4,t 2 = £,4,,tx = £4,43/ = £4, x 2 = £,4,xy = t,4,y 2 = 0, 

£,3,x = £2, y, £3,z + £4, a = 0, £4,1, = £4,11, £2, z + £4,1 = 0, 

and 

Til = -«£l,t + £2,t + <2,y - w£4,^, 
»J2 = ~W^4,y + £3,t - «£2,s, - «£l,i, 

??3 = «£4,x +£4,t + v(,A,y ~ «>£l,t, (4) 

»74 = qrfc/p + 2q^ t , 

V5,t = V5,x = r?5, H = V5,z = r?5,u = V5,v = »75,«i = »75,g = 0, 
T)5 = pri5, P - 

First, equations (3) require that £iS are just 
£1 = ait + a,2, 

£2 = a 3 t + a t x + a 5 y + a 6 z + a 7 , 

£3 = a$t — a^x + ci4j/ + a 9Z + aio, (5) 
£4 = ant — ciax — agy + cuz + ai2, 

where ai,---,oi 2 are arbitrary constants. Then, equations (4) and (5) 
requires that 

771 = (04 - ai)u + a 5 t> — a 6 w + a 3 , 

r/2 = — a&u + («4 — ai)v + (14W + as, 

r]3 = —a^u — agv + (04 — ai)w + an, (6) 

f]4 = (2oi - 2a 4 + 013)0, 

rfs = ai3P, 

where ai3 is an arbitrary constant. Therefore 



3 



Theorem 1. The Lie algebra of infinitesimal symmetries of the system 
(1) is spanned by the 13 vector fields 



(shifts of origin) 



x l 


= d x , 


x 2 


= 9y, 


x 3 


= d z , 


Xi 


= d t , 



Xs = t dx + d u , 

X§ — td y + d v , (uniformly frame motions) 

X 7 = td z + d w , 

(7) 

X s = t d t - u d u - v d v - w d w + 2q d q , 

X 9 = xd x + yd y + zd z + ud u + vd v + w d w — 2r d r , 



Xio = y d x - x d y + v d u - u d v , 

Xn = — zdy +yd z — w d v +vd w , (rotations of reference frame) 
Xi2 = — zd x + x d z — w d u + ud w , 



Xi 3 = qd q + pd p , 

These infinitesimal symmetries will generate a Lie algebra Q over the field 
of real or complex numbers. The commutator table of Lie algebra g for (1) 
is given below, where the entry in the ith row and jth column is defined 
as [Xi,Xj] = Xi.Xj - Xj.Xi, i,j = 1, ■■■,13. 





Xi 


x 2 


x 3 


Xi 


x 5 


x 6 


x 7 


Xs 


x 9 


Xio 


Xn 


Xi 2 


Xl 3 


x 1 




























-x 2 





x 3 





x 2 


























x 2 


Xi 


X 3 








x 3 


























x 3 





-x 2 


-Xi 





x 4 














X! 


x 2 


X 3 


Xi 

















x 5 






















-X 5 X 5 


-x 6 





x 7 





x 6 











-x 2 











-x 6 


X G 


X 5 


x 7 








X 7 











-X 3 
















-X 6 


-X 5 





X s 











-Xi 


x 5 


X 6 


x 7 




















x 9 


-x 1 


-x 2 


-x 3 





-x 5 


-X 6 


-x 7 




















Xio 


x 2 










X 6 


-X 5 














— Xl2 


Xn 





Xu 





-X 3 


x 2 








-x 7 


x 6 








Xi 2 





—Xio 





X\2 


-x 3 





Xi 





-x 7 





x 5 








-Xn 


Xio 








X13 










































3 Reduction of the system 

Now, we look at the group- invariant solutions to the system (1). 

The system (1) is a sub-manifold M of jet space J 2 (R 4 ,R 5 ) defined 
by the following equations 

u t + uu x + vuy + wu z + p x /q = 0, 
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vt + uv x + vv y + wv z +p y /q = 0, 

w t + uw x + vw y + ww z + p z /q = 0, (8) 
qt+q (u x + v y + w z ) + uq x + vq y + wq z = 0, 
p t + jp (u x + v y + w z ) + up x + vp y + wp z = 0. 

Doing as section 3.1 of [6] and find (in a sense) the most general group- 
invariant solutions to the system (1). 

Theorem 2. The one-parameter groups gi(s) : M — > M generated by 
the Xi, i = 1, • • • , 13 are given in the following table: 

91 ■ (t,x + s,y,z,u,v,w,p,q) 

92 ■ (t,x,y + s,z,u,v,w,p,q) 
9s ■ (t,x,y,z + s,u,v,w,p,q) 
94 ■ (t + s,x,y,z,u,v,w,p,q) 

g 5 : (t,x + st,y,z,u + s,v,w,p,q) 
g 6 : (t,x,y + st,z,u,v + s,w,p,q) 

97 ■ (t,x,y,z + st,u,v,w + s,p,q) (9) 
gs : (e s t, x, y, z, e~ s u, e~ s v, e~ s w,p, e 2s q) 
g$ : (t,e s x, e s y,e s z, e s u,e s v, e s w,p,e 2s q) 
gw : (t, x coss + ysins, y cos s — a; sin s, z, 

«sins + u cos s,v coss — u sin s, w,p, q) 
gn : (t, x, y cos s — z sin s, z cos s + y sin s, 

u, v cos s — w sin s, w cos s + v sin s, p, q) 
gi2 : (t, x coss — z sin s, y, x sin s + z cos s, 

wcos s — w sins, w, w cos s + usin s,p, g) 
3i3 : (t,x,y,z,u,v,w,e s p,e s q) 

Where, entries give the transformed point 

exp(sXi)(t,x,y,z,u,v,w,p,q) = (t,x,y,z,u,v, w,p, q). 

Since each group gi, i = 1, • • • , 13, is a symmetry group of (1), then 

Theorem 3. Let i = 1, • • • , 13 ; s € K, and Ael - {0}. // 

u = U{t,x,y,z), v = V(t,x,y,z), w = W(t,x,y,z), 
p = P(t, x,y,z), q — Q(t, x, y, z), 
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is a solution of the system (1), so are the functions 

u (i) = u (i) {t,x, yi z), v (i) = V (l) (t,x,y,z), w w = W (l) (t,x,y,z), 
pW=pW{t,x,y,z), q (i) =QV(t,x,y,z), 



LUflCl fc. 






u« = U(t,x + s,y,z), 


= 


U(t,x,y + s,z), 


v (1) = V(t,x + s,y,z), 




V(t,x,y + s,z), 


W W = W(t,x + s,y,z), 


w m = 


W(t,x,y + s,z), 


= P(t,x + s,y,z), 


p (2) = 


P(t,x,y + s,z), 


— Oii t -1- q ii y\ 

q — d, -f b, y, zj . 


„(2) _ 
1 — 




u (3) = U(t,x,y,z + s), 


ti< 4 > = 


U(* + s,a;,y,z), 


= v(t,x,y,z + s), 


V W = 




w {3 > = W(t,x,y,z + s), 


(A \ 




p (3) = P(t,x,y,z + s), 


p (4) = 




g (3) = Q(t,x,y,z + s). 


,(4) = 


Q(t + s,x,y,z). 


u (5) = U(t,x + st,y,z) - s, 


u< 6 > = 


U(t,x,y + st,z), 


« (5) = V(t,x + st,y,z), 


„(6) = 


V(t, x, y + st, 2) — 


w (5) = W(t,x + st,y,z), 


w (6) = 


: W(t,x,y + st,z), 


p< 5 > = P(t,x + st,y,z), 


p (6) = 


P(t,x,y + st,z), 


q (5) = Q(t,x + st,y,z). 


g (6) = 


Q(t,x,y + st,z). 


u m = U(t,x,y,z + st), 




= X.U(X.t,x,y,z), 


= V(t,x,y,z + st), 


„(8) 


= X.V(X.t,x,y,z), 


w m = W(t,x,y,z + st) - s, 


w w 


= X.W(X.t,x,y,z) 


p {7) = P(t,x,y,z + st), 


p (8) 


= P^X.t, x, y, 2) , 


q (7) = Q(t,x,y,z + st). 


? (8) 


= jz-Q(X-t,x,y,z) 



tt (9) = ±.U(t,X.x,X.y,X.z), u (13) = U(t,x,y,z), 

v (9) = ±.v(t,A.a;,A.j/,A.z), « (13) = V (t, x, y, z) , 

w 0) = j-.w(t,X.x,\.y,\.z), w (13) = W(i,a;,y,«), 

p (9) = P(t,X.x,X.y,X.z), ' p (13) = X.P(t,x,y,z), 

= X 2 .Q(t,X.x,X.y,X.z). g< 13 > = X.Q(t, x,y, 2) . 



u ' — cos s.U[t, x cos s + y sins, y cos s — resins, zj 
— sin s.y (t, x cos s + y sin s, j/ cos s — a: sin s, 2^ 

i/ 10 ' = sin s.ll(t, x cos s + j/sin s, j/coss — :rsin s, 2) 
+ cos s.V(t, a; cos s + 2ysins,ycoss — a; sins, 2 



G 



,(10) 



,(10) 



.(10) 



W(t, x cos s + y sin s, y cos s — x sin s, 2) , 
P(t, x cos s + y sin s, j/ cos s — x sin s, 2) , 
Q(t, 2: cos s + j/ sin s, j/coss — £ sin s, 2) . 



(ii) 
(ii) 



„(") 



(11) 
.(11) 



U(t, x, j/cos s — z sins, 2 cos s + j/sin s), 
cos s.V(t, x, j/cos s — 2 sins, 2 cos s + j/sins) 
+ sin s.lWt, x, j/cos s — 2 sins, z cos s + j/sin s) , 
cos s.W(i, x, j/cos s — 2 sins, z cos s + j/sin s) 
— sin s.V (t, x, y cos s — 2 sin s, 2 cos s + j/ sin s) , 
P(t, x, y cos s — 2 sin s, 2 cos s + j/sin s) , 
Q(t,x,ycoss — 2 sins, 2 cos s + j/sins). 



,(12) 



,(12) 
,,(12) 



,(12) 
.(12) 



cos s.U (t, x cos s — 2 sin s, y, x sin s + 2 cos s) 
+ sin s.lWt, x cos s — 2 sin s, j/, x sin s + 2 cos s), 
V(i,xcoss — 2 sin s, j/, x sin s + 2 cos s), 
coss.lWt, xcoss — 2 sins, j/,xsins + 2 cos s) 
— sin s.U (t, x cos s — 2 sin s, j/, x sin s + 2 cos s) , 
P(i,xcoss — 2 sins, j/,xsins + 2COSs), 
Q(i, xcoss — 2 sins, j/,xsins + 2coss). 



4 Structure of Lie algebra g 

In this section, we determine the structure of full symmetry algebra Q. of 
system (1). 

The center 3 of g is Span R {Xi3}. Therefore, the quotient algebra 
fli ~ 2/3 is Span R {Yi, ■ ■ ■ , Y12}; where Y t := Xi + 3, i = 1, ■ ■ ■ , 12. 
The commutator table of Lie algebra Qi is given below, where the entry 
in the ith row and jth column is denned as = Yi.Yj — Yj.Yi, 



7 



i,j = !,-■■, 12. 





Fi 


Y 2 


Y 3 


Y 4 


Y 5 


Y 6 


Y 7 


Y s Y 9 


Y10 


Yn 


Yl 2 

























Yi 


-y 2 





Y 3 


y 2 























Y 2 


Yi 


Ys 





y 3 























y 3 





-Y 2 


-Yi 


y 4 














Yi 


Y 2 


Y 3 


y 4 











y 5 











-Yi 











-Y B Ys 


-Ye 





Y 7 


y 6 











-Y 2 











-y 6 y 6 


Ys 


Y 7 





Y 7 











-Ys 











-Yr Y 7 





-Y 6 


-Ys 


Y s 











-Yi 


Y 5 


Y 6 


Y 7 














Y 9 


-Yi 


-Y 2 


-Y 3 





-Y, 


-Ye 


-Y 7 














Y 10 


Y 2 


-Yl 








Y 6 


-Ys 











-Yi 2 


Y11 


Fn 





-Y 3 


Y 2 








-Y 7 


Y e 





Y12 





-Y10 


Y12 


-Y 3 





Yi 





-Y 7 





Y 5 





-Y11 


Y10 






The algebra gi is not solvable, because 
fli 1 ' :=[fli,fli] = Span R {Yi ) ...,Yio}, and < 2) := «] = fl« . 



Lie algebra Qi admits a Levi-decomposition into the semi-direct sum 

where r = Span R {Yi, • ■ ■ , Y9} is the radical of Qi (the largest solvable 
ideal contained in Qi) and S = Span R {Yio, Yn , Y12} is a semi-simple 
subalgebra of Qi. The radical r is uniquely defined but the semi-simple 
subalgebra S is not. S as an algebra is isomorphic to 3— dimensional special 
orthogonal algebra so(3) := {A € Mat (3 x 3;R) | A* = -A}, which is a 
simple algebra. 

The radical t is solvable, with and 

r (1) =Span R {Yi,...,Y 7 }, t (2) = Span R {Yi, Y 2 , Y 3 }, r (3) = {0}. 

It is semi-direct sum c = K Q 2 , where Q 2 = Span R {Yg, Yg} is isomor- 
phic to the two dimensional Abelian Lie algebra R 2 ; i.e. 

r = r (1) k R 2 . 

r' 1 ' is semi-direct sum = r (2) IX 03, where 03 = Span R {Y4, • • ■ , Y7} 
is isomorphic to the four dimensional Abelian Lie algebra R 4 , and r' 2 ^ is 
isomorphic to the three dimensional Abelian Lie algebra R 3 ; i.e. 

r (1) =R 3 k M 4 . 

Theorem 4. The symmetry algebra q of system (1) have the following 
structure 

g = R x (((M 3 k R 4 ) K R 2 ) k so(3)). 
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Theorem 5. There is a normal Lie-subalgebra of symmetry Lie-group 
of system (1) which is isomorphic to the Galilean group o/M 4 : 



Gal(4,l 




R G 0(3, M), s G R, and y, v G R 3 



Therefore, the system (1) is invariant up to Galilean motions of space- 
times R 4 (see [3J). 

Proof. Let b = Span R {Xi, ■ • • , X7, X10, Xn, X12}. Lie algebra structure 
of b is as ga[(3), the lie algebra of Galilean group Gal(4, R), and b an ideal 
of g. Therefore, there is a Lie-subgroup of G such that its Lie-algebra is 
b, by the Theorem 2.53 of Olver [6]. □ 

Conclusion 1. q is semi-direct sum of 

ga[(4) ^ Span R {Xi,---,X 7 ,Xi ,Xii,Xi 2 } 
and 3— dimensional Abelian algebra R 3 = SpanjXs, X9, X13}. 



5 Optimal system of sub-algebras 

As is well known, the Lie group theoretic method plays an important role 
in finding exact solutions and performing symmetry reductions of differen- 
tial equations. Since any linear combination of infinitesimal generators is 
also an infinitesimal generator, there are always infinitely many different 
symmetry subgroups for the differential equation. So, a mean of deter- 
mining which subgroups would give essentially different types of solutions 
is necessary and significant for a complete understanding of the invariant 
solutions. As any transformation in the full symmetry group maps a solu- 
tion to another solution, it is sufficient to find invariant solutions which are 
not related by transformations in the full symmetry group, this has led to 
the concept of an optimal system [7]. The problem of finding an optimal 
system of subgroups is equivalent to that of finding an optimal system of 
subalgebras. For one-dimensional subalgebras, this classification problem 
is essentially the same as the problem of classifying the orbits of the ad- 
joint representation. This problem is attacked by the naive approach of 
taking a general element in the Lie algebra and subjecting it to various 
adjoint transformations so as to simplify it as much as possible. The idea 
of using the adjoint representation to classify group-invariant solutions 
was due to [7] and [5] . 

The adjoint action is given by the Lie series 

Ad(exp(sX i )X j ) = Xj - spCi, Xj] + ^[X ; , [X;, Xj]] - ■ • • , (10) 

where [Xi,Xj] is the commutator for the Lie algebra, s is a parameter, 
and i,j = 1, • • • , 13. We can write the adjoint action for the Lie algebra 
Q, and show that 
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Theorem 6. A one-dimensional optimal system of (1) is given by 

1) X = X 8 + ^]i g a i X i , 

2) X = a 7 X 7 + Xs + X 9 + ^] 3 =10 a i X i , 

3) X = a 6 X 6 + arXr + X 8 +X 9 + J2]i 11 aiX i , 

4) X = a 5 X 5 + X 8 + X 9 + auXi + ai 3 X 13 , 

5) X = Y,l= 5 a ^+ X ^+ X a + a ^is, 

6) X = a 4 X 4 + X 9 + J2ll 10 OiXi, 

7) X = X 4 +a5X B + o 8 X 8 + X;" 10 OiX i , 

8) X — Xa + aeXtj + a$X$ + aioXio + arzXrz + 013X13, 

9) X — X4 + 07X7 + a$X 8 + aioXio + 012X12 + 013X13, 

10) X = X 4 + ELb a ^ + a 13Xl3, 

11) x = x 4 + a 7 x 7 + E" 10 "iXi, 

12) X = X 7 + ^ 10 aiX i , 

13) X = a3X 3 + J2ll 10 aiXi, 

U) X = X 4 + a 6 X 6 + J2Zu a > X i> 

15) X — X4 + 05X5 + 011X11 + 013X13, 

16) X = X 4 + Y,L 5 a * x * + 013X13, 
J 7; X = X 6 + X)" u OiXi, 

18) X = a 2 X 2 + ai Xi, 

19) X = a 2 X 2 + 012X12 + 013X13, 

20) X = X s + 011X11 + 013X13, 

21) X = a 2 X 2 + 03X3 + X 5 + 013X13, 
22; X = aiXi + Xn + 013X13, 

23) X = aiXi + 03X3 + X 6 + 013X13, 

24) X = aiXi + 02X2 + 03X3 + 013X13, 

25) X = Xg + aioXio + anXn + 013X13, 

26) X = a 3 X 3 + X10 + 011X11 + 013X13. 

Proof: F(sXi) : Q — > g defined byln Ad(exp(sXi)X) is a linear map, 
for i = 1, • • • , 13. The matrix of F(sXi), i = 1, ■ • • , 13, with respect 

to basis {Xi, • • • , X13} is 

1 0000000000" 

1 000 00 0000 

1 0000000000 

1000000000 

0100000000 

0010000000 

0001000000 

0000100000 

-s 0000010000 

s 0000001000 

000 000100 

O-sOOO 000010 

0000000001 
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and /13, respectively. Let A = cnXi, then 

F(s7.X 7 )oF(s 6 .X 6 )o-..oF(si.Xi) : X i-> 

(s3ai2 — S2ttio — siag + ai — S5S4CI8 + sr,a4 — S4as).Ai 

+ ( — S6S408 + S3«ll — S2CI9 + SiaiO +132+ Sea4 — S406).X2 

+(— S7S4O8 — S3i9 — s 2 au — s 1 a 12 + a3 + s 7 a 4 — s 4 a 7 ).X 3 (11) 
+ (a 4 - s 4 a s ).A 4 + (-s 6 a 10 - s 5 a g + s 5 a 8 + a 5 + s 7 ai 2 ).X 5 
+ ( — sea9 + S6a 8 + ssaio + a 6 + S7an).X 6 

+ (S708 — «6dll — Sodl2 + 0,7 — Sjat^.Xj 

+as-Xs + ag.Xg + aio-Aio + fln.Xn + ai 2 .Xr2 + 013. X13 

If as, and 09 7^ 0, and a 8 7^ 2a9, then we can make the coefficients of 
Xi, - ■ ■ ,Xj vanish. Scaling X if necessary, we can assume that as = 1. 
And X is reduced to Case 1. 

If as — ag / 0, and a w 7^ 0, then we can make the coefficients of 
Xi ■ ■ ■ ,Xe vanish. Scaling X if necessary, we can assume that os = 1. 
And X is reduced to Case 2. 

If as = ag / 0, aio = 0, and a\ 2 0, then we can make the coefficients 
of X\ , • • • , X5, and Aio vanish. Scaling X if necessary, we can assume that 
a8 = l- And X is reduced to Case 3. 

If as = ag ^ 0, aio = a\ 2 = 0, and an 0, then we can make the 
coefficients of Ai,---,A4, Xe, X7, Aio, and X12 vanish. Scaling X if 
necessary, we can assume that as = 1. And X is reduced to Case 4. 

If as = ag 7^ 0, aio = an = «12 = 0, then we can make the coefficients 
of Xi, ■ ■ ■ , X4, and X10, • ■ ■ , An vanish. Scaling X if necessary, we can 
assume that as = 1. And A is reduced to Case 5. 

If as =0, and ag 7^ 0, then we can make the coefficients of Xi, ■ ■ ■ , X3, 
and X5, ■ ■ ■ , Xs vanish. Scaling A if necessary, we can assume that ag — 1. 
And X is reduced to Case 6. 

If as 0, ag = 0, and an 7^ 0, then we can make the coefficients of 
Xi, • ■ ■ , A3, and Xq, Xi, and Xg vanish. And A is reduced to Case 7. 

If as 7^ 0, 09 = an = 0, and ai2 7^ 0, then we can make the coefficients 
of Ai, • • • , A3, and Xs, X7, Xg, and An vanish. And X is reduced to 
Case 8. 

If as 7^ 0, 09 = an = ai2 = 0, and aio 7^ 0, then we can make the 
coefficients of Xi, ■ ■ ■ , A3, and X$, X§, Xg, and An vanish. And X is 
reduced to Case 9. 
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If as / 0, and ag = aio = an = ai2 = 0, then we can make the 
coefficients of Ai, ■ ■ ■ , A3, and Xg, ■ ■ ■ , A12 vanish. And X is reduced to 
Case 10. 

If as = ag = 0, 04 =fc 0, and 010 7^ 0, then we can make the coefficients 
of X\, ■ ■ ■ , A3, A5, Xq, X$, and Xg vanish. Scaling X if necessary, we can 
assume that 04 = 1. And X is reduced to Case 11. 

If 04 = as = 09 = 0, aio 7^ 0, ai2 7^ 0, and a6ai2 + 07010 7^ asan, then 
we can make the coefficients of Ai, ■ ■ ■ ,X&, X$, and Xg vanish. Scaling 
X if necessary, we can assume that 07 = 1. And X is reduced to Case 12. 

If 04 = as = 09 = 0, aio 7^ 0, 012 7^ 0, and 00012 + 07010 = 05011, then 
we can make the coefficients of X\, X2, and X4, ■ ■ ■ ,Xg vanish. And X is 
reduced to Case 13. 

If as = ag = aio = 0, and 04 7^ 0, then we can make the coefficients of 
Xi, • • • , A3, A5, and Xf, ■ ■ ■ , Aio vanish. Scaling X if necessary, we can 
assume that 04 = 1. And X is reduced to Case 14. 

If as = ■ ■ ■ = aio = 0, 012 = 0, 04 / 0, and an / 0, then we can 
make the coefficients of X\, ■ ■ ■ , A3, and As, • • • , Aio vanish. Scaling A if 
necessary, we can assume that 04 = 1. And A is reduced to Case 15. 

If as = • ■ • = 012 = 0, and 04 7^ 0, then we can make the coefficients 
of Ai, • • • , A3, and As, • • • , Aio vanish. Scaling A if necessary, we can 
assume that 04 = 1. And A is reduced to Case 16. 

If 04 = as = 09 = aio = 0, 012 7^ 0, and a 5 an / 06012, then we can 
make the coefficients of Ai, • ■ ■ , A5, and A7, • • • , Aio vanish. Scaling A if 
necessary, we can assume that a6 = 1. And X is reduced to Case 17. 

If 04 = as = 09 = aio = 0, 012 7^ 0, asan = a6ai2, and 06 7^ 0, then 
we can make the coefficients of Ai, • • • , A5, and A7, • ■ ■ , Aio vanish. And 
X is reduced to Case 18. 

If 04 = oo = as = 09 = 010 = 0, 012 7^ 0, asan = 060i2, and ae 7^ 0, 
then we can make the coefficients of Ai, • ■ • , A5, and A7, • ■ ■ , Aio vanish. 
And X is reduced to Case 18. 

If 04 = a6 = as = • • • = on = 0, 012 7^ 0, 05011 = 060i2, 06 7^ 0, and 
a 5 / 0, then we can make the coefficients of Ai, and A3, • • ■ , An vanish. 
And X is reduced to Case 19. 

If 04 = OS = 09 = 010 = 012 = 0, 05 / 0, and an / 0, then we can 
make the coefficients of Ai, • • • , 04, Xe, ■ ■ ■ , Aio, and ai2 vanish. Scaling 
A if necessary, we can assume that 05 = 1. And A is reduced to Case 20. 

If 04 = as = • • • = 012 = 0, and 05 7^ 0, then we can make the 
coefficients of Ai, 04, and As, • • • , A12 vanish. Scaling A if necessary, we 
can assume that 05 = 1. And A is reduced to Case 21. 

If 04 = a 5 = as = 09 = aio = 012 = 0, and an / 0, then we can make 
the coefficients of A2, • • • , Aio and 012 vanish. Scaling A if necessary, we 
can assume that an = 1. And A is reduced to Case 22. 

If 04 = 05 = OS = • • • = 012 = 0, and ag 7^ 0, then we can make the 
coefficients of A2, 05, and A7, • • • , A12 vanish. Scaling A if necessary, we 
can assume that a& = 1. And A is reduced to Case 23. 

If 04 = 05 = oo = 0, and as = • • • = 012 = 0, then we can make the 
coefficients of A4, • ■ ■ , A12 vanish. And X is reduced to Case 24. 
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If an = as = ag = ai2 = 0, aio 7^ 0, and araio 7^ asan, then we 
can make the coefficients of X\, ■ ■ ■ , X4, Xe, ■ ■ ■ ,Xg, and a\2 — vanish. 
Scaling X if necessary, we can assume that as = 1. And X is reduced to 
Case 25. 

If C14 = as = ag = ai2 = 0, 07010 = asan, and aio 7^ 0, then we 
can make the coefficients of Xi, ■ ■ ■ , -X4, X@, ■ ■ ■ ,Xg, and 012 = vanish. 
Scaling X if necessary, we can assume that aio = 1- And X is reduced to 
Case 26. And the Theorem follows. □ 

According to our optimal system of one-dimensional subalgebras of the 
full symmetry algebra g, we need only find group-invariant solutions for 
26 one-parameter subgroups generated by X as Theorem 6. For example, 
as a direct consequence of Theorem 3 and case 1 of Theorem 6, we deduce 
the following conclusion: 

Conclusion 2 Let s, ci, • • • ,ce £ K with a, C2, and c$ =fc 0, and s be 

sufficiently small. Then, if 

u = U(t,x,y,z), v = V(t,x,y,z), w = W(t,x,y,z), 
p = P(t, x,y,z), q = Q(t, x, y, 2), 

is a solution of the system (1), so are the functions 

u — CiCjp ( cos(c3s) cos(c5s) — sin(c3s) sin(c4s) sin(css)) -U(t, x, y, z) 

— c i c 2~ S ( sin(c3s) cos(c 5 s) + cos(c3s) sin(c4s) cos(css)) -V(t, x, y, z) 

+CiC2" s cos(c4s) sm(css).W(t, x, y, z), 

v — CiCj 3 sin(c3«) cos(c4s).?7(t, x, y,z) 
+C1C2 s cos(c3«) cos(c4s).V(t, x, y, z) 
+c{c2 a . sin(c4s).W(t, x, y, z), 

w = —c\c^ s ( cos(c3s) sin(c 5 s) + sin(c3s) sin(c4s) cos(css)) -U(t, x, y, z) 
+ c i c 2 S ( sin(c3s) sin(c 5 s) — cos(c3s) sin(c4s) cos(css)) .V (t, x, y, z) 
+CiC2" S cos(c4s) cos(c5s).W(t, x, y, z), 

P = eg 8 ' .P(t,x,y,z), 

q = Ci c 2 c 6 .Q(t,x,y,z), 

where t = cf .t, and 

x = cl ( sin(css) sin(c4s) sin(c3«) + cos(css) cos(c3s)) .x 
+C2 ( sin(c3s) cos(css) — cos(c3s) sin(c4«) sin(c 5 s)) .y 
— c 3 2 cos(c4s) sin(c5s).2, 

y = — C2 sin(c3s) cos(c4s).a: + (C2S) cos(c3s) cos(c4s).y — C2 sin(c4s).2;, 

z = C2 ( cos(c 3 s) sin(c 5 s) — sin(c 3 s) sin(c 4 s) cos(c 5 s)).:r 
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+c 2 ( cos(c3s) sin(c4s) cos(css) + sin(c3s) sin(css)) .y 

+C2 C0S(C4S) COs(C5fi).2. 



The construction of the group- invariant solutions for each of the one- 
dimensional subgroups in the optimal system proceeds in the same fashion. 
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